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The quantum coherence of electrons can be probed by studying weak local- 
ization corrections to the conductivity. Interaction effects lead to dephasing, 
with electron- electron interactions being the important intrinsic mechanism. 
A controversy exists whether or not the dephasing rate, as measured in a 
weak localization experiment, vanishes at low temperatures. We review the 
non-perturbative analysis of this question and some of the arguments which 
have been raised against it. The compact form of the presentation should 
make the derivation more transparent and accessible for discussions. We 
also compare with recent experiments. 

PACS: 73.23.-b, 72.15.-v, 72.70. +m. 



1. INTRODUCTION 

Quantum coherence and interference are basic principles of quantum 
mechanics as well as essential for the concepts of quantum state engineering. 
Loss of phase coherence results in a suppression of interference and signals 
a transition from quantum to classical behavior. Thus, it is crucial to un- 
derstand how quantum coherence is destroyed in various physical systems. 

Isolated quantum degrees of freedom preserve phase coherence. Inelastic 
interactions with other - quantum or classical - degrees of freedom destroy 
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coherence. For electrons in metals at low temperatures the dominant intrin- 
sic interaction is electron-electron scattering. As the temperature is lowered 
further, these inelastic processes freeze out, and the scattering rate vanishesil, 
1/rin oc T^. Largely because of this property it is frequently assumed that 
also the electron decoherence rate I/tj^ should vanish as the temperature 
approaches zero. 

The electron decoherence can be probed by studying the weak localiza- 
tion correction to the conductivity. While many of the earlier experimental 
appeared to confirm the theoretical expectationdlfl'i, more recent workiS 
(see Ref. ^ for further references) suggested that the electron dephasing time 
T^p in disordered meljals saturates at sufficiently low temperatures. Extensive 
experimental checksQ allowed to rule out various sources for this saturation, 
such as electron scattering from magnetic impurities, heating and exter- 
nal high frequency noise. Several other saturation mechanisms had been 
explored, such as 1// noise, interaction of electrons with nuclear magnetic 
moments or with two-level systems in the two-channel Kondo regime!. They 
appear not to play a role in the experiments described in Ref. |^ nor in the 
more recent studies of Ref. ^. As concluded in Ref. |7|, these recent experi- 
mental findings "reinforce the earlier conclusion that the saturation is a real 
effect, most likely arising from electron-electron interaction". 

A low-temperature saturation of t<^ due to electron-electron interactions 
has drastic physical consequencesBily. One of them is that the Coulomb in- 
teraction precludes (strong) Anderson localization in disordered conductors. 
Another is that the Fermi liquid theory may not apply to disordered con- 
ductors. Dephasing at low temperatures also imposes practical limitations 
on quantum state engineering with electronic devices. All that makes the 
low-temperature electron decoherence one of the most intriguing issues of 
modern condensed matter physics. 

In this context it is important to note that several model systems where 
a quantum particle is coupled to a bath (e.g. a bath of harmonic oscillators 
as in the Caldeira-Leggett model§) yield results wh ich indicate a nonva- 
nishing dephasing rate even at zero temperaturcll3ll3ll-il't!j. We further like 
to mention the work of A.G. Aronov and P. Wolflelll which emphasizes the 
fundamental importance of dephasing due to fluctuating magnetic field for 
such issues as high temperature superconductivity and the metal-insulator 
phase transition. 

The aim of the present paper is to review the theoretical analysis of the 
Coulomb- interaction-induced quantum dephasing of electrons in disordered 
conductors derived before by two of the present authors (GZ)t3't3ll2. This 
work has stimulated a heated discussion and exchange of lengthy publica- 
tions. One of our goals is to present in a compact form the main steps of 
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the derivation, which should help an interested reader to judge the validity 
of this approach. Some o f th e more recent work - in particular we like to 
mention that of von Delftt^l - have put the discussion on a more rational 
basis by explicitly pointing out approximations which had been used. The 
second goal of this article is to justify these approximations. Finally, we will 
compare with some recent experiments with the aim to judge whether they 
can help to discriminate between different theoretical views. 

The structure of the paper is as follows: We will first review the deriva- 
tion of weak localization effects on the conductivity as well as a simple anal- 
ysis of the influence of electron-electron interactions. At this level we repro- 
duce some important results, but, in addition, we introduce the formalism 
and notations which will be used in later sections, where the same questions 
will be analyzed on a more rigorous level. We then summarize in cor n pact 
form the derivation which has been presented before by two of us (GZ)li3'llZHl3. 
This part should make explicit which approximations have been used. We 
then will present several arguments which had been raised against the con- 
clusions of GZ as well as our counterarguments. The final section is devoted 
to a comparison with recent experiments. 



2. WEAK LOCALIZATION BASICS 

In this section we shortly review the derivation of the weak localization 
correction to the conductivity for non-interacting electrons using a path 
integral formulation, and we present a hand- waving approach to account 
for the effects of electron-electron interactions. 

The density matrix of a quantum particle evolves according to 

p(tf,rif,r2f) = J (irii(ir2i J(tf,ii,rif,r2f,rii,r2i)/9(ti,rii,r2i), (1) 

where the time evolution operator J combines the forward and backward 
propagators typical for the density matrix. For a single particle which moves 
in an impurity potential U (r) and which is subject to a static applied electric 
potential Vx, this operator can be expressed as a path integral over two paths 

J(tf,ii,rif,r2f,rii,r2i) = 

ri(tf)=rif r2{tf)=r2f 

Vri / Dr2 e •'^^ , (2) 



ri(ii)=rii r2(ti)=r2i 



where So[r] = dt' [mt^{t')/2 - U{r{t'))]. 
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We are interested in the conductivity in spatially homogeneous and 
stationary situations. We, therefore, assume Vx{r) = — Er and calculate the 
current at r and tf due to the field at r(t') and t' in an expansion of the 
density matrix to linear order in Vx- Thus we get 



rif r2f 



a 



3m 



J dt' J driidr2i[Vrif - Vr^f] J Vvi J Vr2 [ri(t') - r2(t')] 



X e 



rii r2 

i5o[ri]-'t5o[r2] 



2J /)(ti,rii,r2i) 



(3) 



rif=r2f=r 



In the stationary limit the result does not depend on time. Accordingly we 
can suppress the dependence of a on the time tf and set t^ —oo. We 
also assume the current to be averaged over the final position rif = r2f = r. 
Hence, the conductivity @ does not depend on this coordinate. 

When the initial density matrix is chosen to be the equilibrium form pQ 
of non-interacting electrons - which does not evolve in time - the expression 
for the conductivity can be simplified further. Since Eq. (|3|) contains no 
time dependence for times earlier than t' , the path integral can be restricted 
to times later than t' . Hence we get 



3m 



J dt' J dr'^dr'2 [Vrif - V,.2t] Jo{t{, t',rif, r2f , r'^, r'; 



X [r'l Po(ri,r'2^ 



(4) 



rif=r2f 



where the time evolution operator Jq is obtained from Eq. (|2|) by setting 
Vx = 0- The expression (Q) can be shown to be equivalent to the standard 
Kubo formula for non-interacting electrons (see e.g. Ref. |20| ). 

In metals, where the electron wavelength is short, one can evaluate 
the path integral Jq in the 'quasiclassical' approximation. This assumes 
that the main contributions arise from classical paths r'^'(t) and quadratic 
fluctuations around them. In the present problem we have to specify the 
paths on both branches of the 'Keldysh contour', rf{t) and rf{t). The main 
contributions to the conductivity (^) arise from pairs of identical classical 
paths, rf{t) = rfit). They yield the Drude conductivity cjd = 2e'^NoD, 
where D is the diffusion constant, proportional to the elastic mean free 
time Tq. As shown by Altshuler, Aronov and Khmelnitsky (AAK)a and 
emphasized by Chakravarty and Schmid (CS)i, the first quantum correction 
to the conductivity, due to 'weak localization', is derived from pairs of time- 
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reversed classical paths rf{t' + t) = rf{t{ — t). They yield the correction 

Sct'wl = -'^ I dtW{t), (5) 
vr J 

Te 

where W{t) = {ATTDt)~'^/'^ a'^~^ is the classical probability for a diffusive 
particle to return to an initial point after time t. Here d is the dimensionality 
of the sample and a the film thickness for d = 2, or the square root of the 
wire cross section (area) for d = 1. The weak localization correction diverges 
in one and two dimensions due to the contributions from the upper limit of 
the integration. This divergence is removed, when interactions are taken into 
account. Their effect is to cut the integral at times exceeding the 'dephasing 
time' T^. This time should be determined in the following. 

Under ideal situations at low temperatures the dephasing time is 
limited by electron-electron interactions. In order to determine it we proceed 
in two stages. We first present a hand- waving approach in the spirit of AAKH 
and CS§, which is accepted to be sufficient at not too low temperatures. In 
the next section we will substantiate and improve the derivation. 

The essence of the simple approach is to replace the interaction by 
an equivalent Nyquist noise source. I.e. the electrons move in the effective 
fluctuating Gaussian potential 6V(t, r) with the correlator 

r,..,„,o), . D . / ^ »"*"*^'^' 

which is fixed by the fluctuation-dissipation theorem. Here e{uj, k) is the 
dielectric^usceptibility of the material. The fluctuating potential should 
be addeda to the action of the electron, 50^5" = /^*' dt — U{r) + 

eSV{t,r)], which enters, e.g., the expression for the conductivity (|3[). 

In the presence of the Nyquist noise a reduction of the time span of the 
path integration, as done in the transition from (^) to (^), is no longer rigor- 
ously allowed since the density matrix now depends on time. Any arbitrary 
density matrix at ti = — oo should relax to the equilibrium density matrix of 
the interacting problem. The use of a reduced form equivalent to (Q) with, 
e.g., the equilibrium density matrix of the non- interacting problem at time 
t' , in general may lead to spurious results describing the decay of this spe- 
cific initial state. However, as will be shown below, for the weak localization 
correction this approximation is justified with sufficient accuracy. 

Due to the Gaussian nature of the fluctuating field 6V one can perform 
the average and finds that the exponent in Eq. (P) acquires an imaginary 
part, 

gi5o[ri]-iS'o[r2]-Si[ri,r2] 
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Fig. 1. A pair of classical paths which contributes to the weak localization 
correction. 

with 

2 *f *f 

Si = ^ JdtJ di [I{t - i, ri(t) - ri(f)) + I{t - i, r2{t) - r2{i)) 
ti ti 

-I{t - i,ri{t) - r2(t)) - I{t - i,r2{t) - ri{i))] . (8) 

The effects of interactions, which are responsible for the decay of cor- 
relations and coherence, are lumped in the imaginary part of the action Si. 
As long as the interactions are weak we can analyze their effect within the 
quasiclassical approximation. For this we need to evaluate Si on the im- 
portant classical paths r2(t) and ri(t), satisfying the classical equation of 
motion without interaction 6S()[r]/5r = 0. Since we are interested in weak 
localization corrections, we concentrate on paths involving pairwise time- 
reversed paths. The leading contribution arises from pairs of paths as shown 
in Fig. 1. The two paths are equal, r2{t) = ri(t), for times earlier than t' , 
they are time-reversed, r2(t) = ri(tf + t' — t), for times t' < t < if, and they 
pass the same region within a diameter of order of the mean free path I at 
two times, t = t' and t = tf. In principle, there could be additional loops 
with time-reversed parts before t' . However, we can ignore such contribu- 
tions because they constitute higher-order weak localization corrections to 
the conductivity, which decay faster in time. For the paths considered, it is 
easy to see that the double time integral in Eq. (^) reduces to 

fk rk ftf rtf 

/ dt dt^ dt dt . (9) 

Jti Jti Jt' Jt' 

The imaginary part of the action 5i(tf — t') grows with the time difference 
tf — t' . For a given pair of paths it may fluctuate on a short time scale of 
order Tg, however this averages out when we sum over all possible paths. 
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The classical paths of electrons in disordered metals are diffusive. When 
we add the contributions of all possible pairs of time-reversed paths a de- 
caying factor {W{t{ — t')e~'^^(*f~* •*) appears under the integral (§). Our aim 
is to find the characteristic time scale r<^ characterizing this decay. In or- 
der to do that, we may replace {W{t)e~^^^^) by W{t)e~^^^^^^K This simple 
averaging has been also performed by CScI, while AAKS evaluated directly 
{W{t)e~^^^^^). However, both approaches gave the same expression for 
up to a prefactor of order one, which is of little importance. 

After these steps the weak localization correction to the conductivity 
(|5|) is modified to become 

5awL = -^^ / dtW{t)e~^''^'\ (10) 
vr J 

Tc 

where fdit) = {Si{t)) has the form 
t t 

fa{t) = e^ J dt' J dt" J d'^xD{\t'-t"\,x)[I{t'-t",x)-I{t'+t"-t,x)] . (11) 



It depends on the 'diffuson' D{t,r) = {AnDty'^/^a'^-''^ exp{-r^ /2t), i.e. the 
probability that an electron propagates the distance r within a time t. 

Using the Drude formula for the dielectric susceptibility of disordered 
metal, e(w, k) = 4iT/{—iLo + Dk'^), we can express the function fd{t) asEl 

3d-d^-2 o^.„^i_d , , , „ / 

U) — UJ 



2 



(cj — oj') coth ■ 



1 — COSCJt UJ COS LOt — COS Uj' V 

LO coth 



u;2 2T uj''^ - u;2 



2T 

(12) 



It is obvious that fd{t) grows with time and limits the range of time integra- 
tion. TheJephasing time can be defined by the condition fd{Tip) ~ 1, with 
the resultty 

1 Ae'^D r^/^- duj f d'^q a; coth(u;/2r) 



cTBa^-'i Ji/r^ 27r J (27r)rf tu^ + (Dq^)^ ' ^^^^ 
Explicitly we find0§ 

=^^[1 + 2T^ in ID, 
l/r^ = [1 + 2rre ln(rr,)] in 2D, 

1/V = .^, %,,A l + 6iTrem in 3D. (14) 
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The divergence of the integral in ([T^) at high frequencies is cut at u; ~ l/xe 
because a classical path needs a time exceeding Tg to return to the same point. 
Because of the quantitative ambiguity related to this cut-off procedure, as 
well as the definition of 1/t<^, the results are valid only up to numerical 
prefactors of order one. 

We observe that the dephasing rate (|l^) does not vanish even at zero 
temperature, and, in this respect, differs from the earlier exnectations. The 
main difference between our work and that of AAKS and CScI lies in the fact 
that they ignore fluctuations of the potential with frequencies higher than 
T. This procedure yields a dephasing rate which vanishes at T = 0. For a 
comparison of the results it is useful to split (^) as 

^-^ + 7(ry ^^^^ 

Here r^o is our - controversial - zero temperature dephasing time, and t{T) 
is basically the result derived by AAKeI. The latter has also been supported 
by Aleiner, Altshuler and Gershenson (AAG)0. Within the framework of 
first order perturbation theory in the interaction they arrive at 

1 _ e^D r duj_ r d'^q u)[coi\i{uj /2T) + tanh((^ - uj)/2T)] 
^ ~ Taj^a^-'i J 2^ J (2^ ^ {u^ + (Dq^)^) cosh^iC/2T) ' ^ ^ 

In this form, after performing the integration over ^ one observes that high 
frequencies w > T do not contribute to the integral, and (^) effectively 
reduces to the result obtained by AAKQ. Below we will comment on the ar- 
guments for either of the results. Before doing so we present, in the following 
section, a more thorough analysis of interaction effects. 



3. INTERACTION EFFECTS: MORE RIGOROUS 

APPROACH 

In this Section we review a more thorough description of interaction 



effects in disordered conductors, based on the work of Refs. 17 and |lq, which 
takes into account Pauli principle etc. The formulation of the problem of 
electron-electron interaction in metals is similar in many respects to the 
problem of a particle interacting with a bath of harmonic oscillators, which 
is known as the Caldeira-Leggett (CL) model. 

First, the fluctuating electro- magnetic field needs to be quantized. This 
is achieved by introducing two independent fluctuating potentials Vi and V2 



for the two branches of the Keldysh contour (see, e.g., Ref. 22) or, equiva- 



lently, the combinations = {Vi + V2)/2 and V = Vi — ¥2- We assume 
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the fluctuations to be small, and expand the action up to the second order 
in V^. This approximation is equivalent to RPA. The first order term in 
the fields describes the motion of electrons in a effective static potential of 
other electrons. It leads to the so-called Hartree diagrams in perturbation 
theory with closed electronic loops. In a homogeneous Jellium model these 
terms vanish because of charge neutrality, but in real disordered metals this 
is not the case. Nevertheless, these terms are not important for dephasing, 
and we ignore them, leaving in the action only the terms quadratic in V^. 
Then the properties of the electro-magnetic field are fully determined by the 
correlators of V^, i.e. the photon Green functions. In homogeneous cases 
they have the form: 

(y+(t,r)F+(0,0)) =/(t,r), {V+{t,r)V~{0,0)) =iR{t,r), (y-y-)^O, 

(17) 

where 

and the function /(i, r) is defined by Eq. (|6|). Note that the function R{t, r) 
vanishes identically for t < due to causality. 

The density matrix of an interacting electron in a metal, />(ri,r2) = 
(^'^(r2)^'(ri)), can be expressed in terms of the density matrix /9y(ri,r2) 
of an electron in the fluctuating potentials as follows: 

p(ri,r2) = (/9y(ri,r2))y± . (19) 



The averaging over V is performed with the aid of (17). The matrix 



pi/(ri,r2) satisfies the generalized nonlinear Liouville equation: 



I- 



[Ho-eV+,pv]--{l-pv)eV~pv--pveV-{l-pv) . (20) 



Here and below in Eq. (21-26) the products are understood in an opera- 
tor sense and integrations over internal spatial (not time) coordinates are 
assumed. Eq. (^) is exact, accounting, for instance, for the Pauli prin- 
ciple. The properties of the interaction are determined by the correlators 
(17). Eq. (pO|) had been derived in Ref. lO, and below we will outline this 



derivation. We can add that Eq. (|2[)D has been independently re-derived by 
Eriksen et alE3 and also by von DelftEJ, who applied Wick's theorem. 

The full 2x2 matrix Keldysh Greens function Gy(ii, t2) rij r2) in the 
fluctuating potentials V = ViSij (for i = 1, 2) satisfles the equation 



{iid/dh - IHo + eV)Gv = aJih - t2)6{ri - ra) . (21) 
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It follows from Eq. (piD that the matrix elements of Gy can be expressed by 



the evolution operators Ui{ti,t2) = Texp —i J^^ dt (Hq — eVi 



as follows: 



Guy = -ie{ti-t2)Ui{tiM)+iUi{tut)pv{t)Ui{t,t2), 

G22,V = -iO{t2-ti)U2{tl,t2)+iU2iti,t)pvit)U2{t,t2), 

Gi2,V = iUi{ti,t)pv{t)U2{t,t2), 

G2iy = -iU2{ti,t)[l- pv{t)]Ui{t,t2). (22) 



On the other hand, G satisfies the Dyson equation 

Gv{ti,t2) = Go{ti,t2) - f' dt'Go{h,t')eV{t')Gv{t',t2), 



(23) 



where the Keldysh matrix Go is defined by Eq. ( |22D with Vi = 0. Eq. (^ 
is compatible with Eq. ( p3[ ) only if 



pvit) = Ui{t,ti){l - Uiit,ti)po[Ui{ti,t) - U2iU,t)]} ^poU2{ti,t), (24) 

where po is an arbitrary time-independent initial density matrix. One can 
verify directly that the density matrix (24) satisfies the equation of motion 



(pOj). This completes the derivation of (|20|). 

The solution ( p^ ) in principle allows one to derive an exact expression 
for the conductivity. Alternatively, we derive here the conductivity directly 
from (pO| ). In the presence of a weak external electric potential, Vx{r) = — Er, 
Eq. (pO|) can be linearized in Vx to yield 



.d6p 



V 



dt 



Hi6pv - 5pvH2 - [eVx,pv], 



(25) 



where 6pv is proportional to Vx and 



Hi = Ho- eV+ - i(l - 2pv)eV-, H2 = Ho - eV+ + \eV~{l - 2py) -(26) 
Continuing along the lines outlined in Section 2, we find the conductivity, 



a 



e 



J dt' J dr[dr2 [Vr-^f - Vr^^] 
-00 

x( Jy(tf,t',rif,r2f,r'i,r2) [r'^ - pvit' ,r[,r2] 



(27) 



rif=r2f 



where 



Jy(t,t',rif,r2f,r;,r'2) = J PpiPriPp2Pr2 e^^^lPi'^i'^^J-^^^IP^'^^'^*], (28) 
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and 

if 

S'1,2 ^ J 



The only approximation used so far is RPA, which ahows a restriction to 
Gaussian fluctuations of the potentials V^. If one averages over with the 



aid of full nonlinear action of the electro-magnetic field, then Eqs. (27-P9|) 
are exact. It is also useful to note at this point that the form of the action 
(|29| ) is not unique since the Hamiltonians (p6|) are sensitive to the ordering of 
momentum and coordinate operators. Under these circumstances details of 
quantization procedure become important. However, all these uncertainties 
are of quantum origin and disappear in the quasiclassical limit ^ ^ 0, which 
will be the only interesting limit in what follows. 

From here on we resort to further approximations. First, in Eq. ( [27| ) we 
replace pv by the equilibrium density matrix of non- interacting electrons pQ . 



As has been pointed out in Ref. 19 this is equivalent to the assumption of a 
factorized initial density matrix in the CL and similar models, an approxi- 
mation which is sufficient in many cases, although for certain correlators it 
is known to fail in the long-time limitlli. Second, also in Eq. ([29| ) we replace 
Pv by pq. It is possible to show that for the problem of weak localization 
both these approximations are justified. We will return to this point in the 
discussion in Section 4. Here we only like to note that this statement can be 
justified by rewriting the exact solution ( p^ in the form of a power series, 
each term of which can be handled without replacing pv by po- 

Having made these approximations, we perform the averaging over 
in Eq. (|^) and get the conductivity: 



2 



3m 

—00 



J dt' j dridr2[Vrif - VrgJ J(tf, t' , r^, r2f, r'^, r's) 



rif=r2f 

X [r'l -r2]/9o(r'i,r2) , (30) 



where 



J(tf,t',rif,r2f,r'i,r2) = j VpiVvi j Vp2Vv2 

^ gi5o[pi,ri]-j5o[p2,r2]-iS'R[pi,ri,p2,r2]-Si[ri,r2] ^^^^ 

and S'o[p,r] = /^*/ dt[pr — p^/2m — U{y)]. The parts of the action iSyi and 
describe the effect of the bath of all electrons on the motion of a single 
electron. We do not give explicit expressions for Sr, since it vanishes for 
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time-reversed diffusive paths (see Ref. ^ for details) and thus is irrelevant 
for dephasingB. Here we would only note that this part of the action contains 
the combination 1 — 2po(p>r) = tanh[(i/o(Pjr)— ;u)/2r], which is responsible 
for the Pauli exclusion principle and generates the term with tanh[(^— w) /2r] 
in first order perturbation theory ([l6|). The imaginary part of the action Si 
is given by Eq. (|8|). It is positive for all possible pairs of paths ri and 
r2, except for two identical ones when it is zero. Only contributes to 
dephasing. In the previous section we have shown that it leads to the decay 
of correlations oc exp(— /^(t)), where fd{t) is given by Eq. ([T^). 

The derivation of Eq. (|^) does not involve disorder averaging, except 
that we have used a translationally invariant form of the voltage correlators 
(|T8|), which is appropriate after averaging. On the other hand, one can 
demonstrate that the result ( |30| ) remains valid in the general case when the 
correlators depend on two arguments separately. We also like to point out 
that the formula (|3^) is very similar to the evolution equation for the density 
matrix of a particle interacting with the bath of oscillators (CL). 



4. DISCUSSION 

The derivation presented above leads to the conclusion that the dephas- 
ing time of electrons in disordered metals, as measured in weak localization 
experiments, remains finite at low temperatures. This result deviates from 
the earlier theoretical understanding and accordingly has been heavily crit- 
icized. In this section we summarize several of the arguments presented 
against this conclusion and our counterarguments. This discussion should 
further clarify under which circumstances the present result applies. 

First-order perturbation theory suggests vanishing T = dephas- 
ing. 

Based on a first-order perturbation expansion AAGcll concluded that 
the dephasing rate is given by expression (|l^). It contains the combination 
coth(u;/2T) + tanh[(^ — oj) /2T], which vanishes in the limit ^, T — > 0. 

In contrast, the result (|T^ ) contains only coth(a;/2T) and remains finite 



for T — > 0. The term tanh[(^ — uj)/2T] in (16) comes from the real part of 
the action 5r,, which vanishes for time-reversed pairs of classical paths, and 
thus is of no importance for the decay of coherence in the long time limit. 

Perturbation theory probes only short times. In general it cannot be 
used for the calculation of t^, since this would require an extrapolation to 
long times, e.g., by exponentiation of the first order perturbative result. 
Such a step is ambiguous, as illustrated by the well-studied problem of the 
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quantum decay of a metastable state in the presence of dissipation. The 
semi-classical decay rate is F = ii£xp(— A). The exponent A = Aq + r]Ai is 
the action on a saddle-point patht^. Here the friction coefficient tj describes 
the interaction with a dissipative bath, and ^0,1 (as well as Bq i below) are 
positive and independent of rj. Fluctuations around the saddle-point path 
contribute to the pre-exponent B. In the weak damping limit one getsll3 
i? ~ i?o + V^i- It is obvious that this result and the specific separation 
into prefactor and exponent - although it can be expanded in - cannot 
be derived unambiguously from a purely perturbative calculation. Similarly, 
for the calculation of the dephasing rate it is important to separate the 
pre-exponent, which is of little importance, from the contributions to the 
exponent. In contrast to the first order calculation the saddle-point analysis 
of gzbEI provides a definite prescription for this separation. 

When comparing saddle-point and first-order calculation one should 
note that the former concentrates on contributions from time-reversed pairs 
of classical paths whereas the latter effectively sums up contributions from 
all paths. Therefore, after the saddle-point approximation is performed one 
cannot recover any more certain details of the first-order calculation. 

We, furthermore, would like to point out that even in the framework 
of perturbation theory, after exponentiation one would find a non-vanishing 
- though small - dephasing rate at T = 0. It is not contained in Eq. ( p!6|) 
which is based on Fermi's golden rule. However, already first order peitui- 
bation theory in the interaction yields additional non-golden-rule termfilEJ. 
These terms arise due to the non-Markovian nature of the interaction, as 
a result of which the usual Golden-rule replacement of [1 — cos{ujt)]/uj'^ by 
a (5-function fails at low temperature. The extra contributions are usually 
small compared to the Golden-rule terms, but they are the leading ones 
when the combination of coth and tanh terms cancel. These extra terms 
also contribute to dephasing at T = 0. 

In the ground state there is no electron scattering. 

This argument is based on the known properties of inelastic scattering. 
As long as the energy of an incident particle is lower than that of the first 
excited state of the scatterer, inelastic scattering is forbidden. The inter- 
nal ground state wave function of the scatterer manifests itself only through 
a form- or Debye- Waller factor. One may conclude then that at low tem- 
peratures all electrons are in the ground state, therefore electron-electron 
scattering is suppressed and no dephasing occurs. 

This argument is not conclusive for several reasons. First, if one would 
describe electron-electron interaction in terms of inelastic scattering one 
should take into account that the scatterer is a macroscopic interacting many 
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particle system with a very small spacing between the energy eigenstates. 
The temperatures where the experiments have been performed are much 
higher than this level spacing. 

Second, since the Coulomb interaction is long-range each electron in- 
teracts with many other ones. A perturbative electron-electron scattering 
picture is actually not appropriate. A more appropriate description is in 
terms of collective electron degrees of freedom, for example plasmons. Also 
the Nyquist noise description used here accounts for collective properties. 
Within this picture all electrons form the environment for a given electron. 

Quasiparticles are fully coherent. 

First we recall that the existence of exact quasiparticles is not proven, 
especially in disordered metal. Even if well-defined quasiparticles do exist, 
it would not rule out dephasing of real electrons. It is the electric current of 
real electrons which is measured. An electron can be understood as a combi- 
nation of an infinite number of coherent quasiparticles (provided the latter 
can be constructed). Therefore it may be only partially coherent even at 
T = 0. A similar scenario can be analyzed rigorously in a Caldeira-Leggett- 
type model, where one (singled-out) oscillator is coupled to an infinite bath 
of oscillators. The single oscillator looses phase coherence even though the 
eigenmodes of the total system are fully coherent. 

During disorder averaging Hikami boxes were missed. 

AAGlill argued that careful disorder averaging within perturbation the- 
ory introduces certain contributions, the 'Hikami boxes'. According to them, 
Hikami boxes correspond to some paths in our approach, which we have 
missed. This would be responsible for the finite dephasing at T = 0. 

On one hand, we understand that Hikami boxes do not correspond to 
any classical path. They appear as a formal feature of the perturbative many- 
body approach used by AAG. There Hikami boxes are needed to restore 
the causality principle. In the path integral approach causality is always 
maintained, and Hikami boxes, together with unusual paths, do not appear 
at all. Moreover, even within the approach of AAG Hikami boxes appear only 
in the terms generated by the action S'r, which are irrelevant for dephasing. 

On the other hand, within our approach the dephasing rate is found 
to be practically insensitive to the particular averaging procedure employed. 
For any given impurity configuration we find that for all time-reversed paths 
the term Si{t) is positive and grows with time. Thus for any pair of such 
paths quantum interference is suppressed by a factor ~ exp(— S'i(t)). Differ- 
ent procedures of averaging over these paths or the impurity configurations, 
e.g. {W{t)e-^^^^^) or I^(t)e-<'^iW> , should yield essentially the same value of 



Low-temperature Dephasing 



dephasing rate, possibly differing by prefactors of order one. 



The replacement pv po is uncontrolled. 

In Ref. it was argued that the replacement pv —>■ Po in Eq. ( p7| ) 
is dangerous. This replacement can be interpreted as preparation of an 
artificial initial state of the system. It might lead to spurious relaxation 
processes, which could be interpreted as dephasing at T = 0. 

When analyzing this question it is important to note that we do not 
study the time evolution of the density matrix of an electron, rather the con- 
ductivity of the electron system. Both evolve differently as can be seen from 
the comparison of Eqs. (||) and (Q). The initial density matrix po(rii,r2i) in 
Eq. (II) is replaced by the combination {r[ — r2)/0o(r'i, in Eq. @). The 
latter can be interpreted as an effective density matrix, which is strongly non- 
equilibrium due to the presence of the factor (r'j^ — r^. This 'non-equilibrium' 
form is the main reason of relaxation and dephasing, rather than the as- 
sumption of factorized initial conditions. More generally, one observes that 
in dissipative systems most of the correlators of the type {A{t)A{0) — ^(0)^) 
decay with time even at T = 0, which does not contradict to the fact that 
the equilibrium density matrix does not relax. 

In spite of these arguments, one might still worry that the replacement 
Pv — > Po could lead to an incorrect long-time asymptotic and wrong conclu- 
sions about the magnitudejif t^. This indeed happens for certain correlators 
in exactly solvable modelai3. We find, however, that it is not the case for 
the conductivity in the weak localization problem. We have already shown 
in Sec. 2 that r^p is not sensitive to the replacement if we can ignore the field 
V~ . The latter can be ignored since it is responsible only for the appearence 
of the part 5r in the action, which gives zero contribution to We can 

further add that even if we include V~ , the replacement pv Po is justi- 
fied. The reason is that the integration in Eq. (27) is confined to the region 
|r'^ — rgl < I. In this range the effect of interactions on the density matrix is 
negligible because I is much shorter than the dephasing length L^p. Finally, 
we have developed an alternative approach based on a rewriting of the exact 
solution (24) in the form of a series. Within this approach the replacements 
Pv — > Po are not made, however the final result for t^^ remains unchanged. 



High frequency modes can be ignored at T = 0. 

The difference between the present and earlier results for the dephasing 
rate lies to a large extend in the inclusion or truncation of high frequency 
modes. Their relative role can be demonstrated if we consider a model 
problem of a particle interacting with a bath of harmonic oscillators. An 
example of such a system is a tunnel junction with high resistance in a 
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circuit involving a dissipative environment. The properties of such a system 
are described by the so-called P(£')-theoryii. The function P{E) is the 
Fourier transform of the self-correlation function of the operator exp^iip), 
with (p = dt'eV{t') being the phase of the junction: 

P{E) = ^ J dt e-^W+^^*, e-^W = (e^^We-^^(°)). (32) 
The function J{t) can be found exactly, J(t) = —F{t) — iK{t), where 

F t) = — / (i6jRe Z(w wcoth ^ , (33 

vr JO 2r iij^ 

K{t) = - dujRe[Z{uj)]-^^. (34) 
vr Jo ^ 



Here Z{lo) is the total impedance of the system tunnel junction plus envi- 
ronment. Note that the function F{t) is similar, although not identical to 
the function fdif) introduced in (^). It is also possible to show that F{t) 
comes from the imaginary part of the action of the environment ^i, while 
K{t) corresponds to the real part S'r. 

We can rewrite the function F{t) approximately as a sum of two terms, 

F(t)^2T— dujRe[Z{uj)] ^ + — / dujRe[Z{uj)] . 

TT Jo UJ^ TT Jo 00 

AAK and other authors considered only the first, 'thermal' term in the 
theory of weak localization. The second term is 'quantum'. It contains the 
contribution of high frequencies. 

We first consider the properties of the function F{t) at T = 0, which 
depend strongly on the impedance Z{ijj). If Re[Z(a;)] = below certain 
frequency Wmin; then exp[— F(t)] first decays with some characteristic time 
Tip, approaches a constant at time t ~ l/wmin, and beyond this time the 
correlations do not decay any more. Now the role of different frequencies 
on the process of dephasing becomes clear. The high frequency properties 
of the impedance determine the initial stage of decay and introduce the 
characteristic decay time r^p, while the low frequency modes determine the 
long-time behavior of the correlations. If the low frequency modes are absent, 
then the decay of correlations stops at time t ~ l/wmin- 

Finite temperature do not change this picture qualitatively. The de- 
cay becomes stronger than at T = 0. At t > l/wmin the decay of corre- 
lations stops again, but the asymptotic value of the correlation function 
(ev5(*)e-*'^(°)) is smaller than for T = 0. 
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On the other hand, if the function Re[Z(a;)] is non-zero at smah uj, the 
correlations decay completely in the long time limit at any temperature in- 
cluding T = 0. 

Further arguments. 



Since the publication of Refs. 10,17 many more arguments have been 



put forward challenging the conclusions reached in these articles. To list a 
few of them: 

- The Pauli principle is not taken into account properly. 

- The electric field is not quantized. 

- The effect of interaction on classical paths is ignored. 

- Detailed balance is violated. 

- The term z5r cancels Si for T ^ 0. 

Here we only want to stress that we have checked these as well as other 
arguments and found in all cases that they did not apply or are in error. 



5. EXPERIMENT 

In our earlier work we have already demonstrated that our theory fits 
the data of Ref. |^. Here we want to add a comparison with some more 
recent experiments. Natelson et alB found that their data for are not 
consistent with the geometry dependence predicted by the standard theoryH, 
and for wide samples they observed saturation of r<^ at low temperatures. In 
the work of the SaclavzMichigan collaboration no saturation was detected in 
silverc3 and pure goldc^ samples down to ~ 50 mK, but a clear saturation was 
observed in copper at temperatures one order of magnitude higher. We will 
show that these results can be explained within the framework of quantum 
dephasing due to electron-electron interactions. 

First we note that Eqs. (14) are valid only in strictly ID, 2D or 3D 



geometries. In practice, however, the diffusion depends on the geometry of 
a particular sample and intermediate regimes may occur. For instance in 
the experiments of Refs. the samples are quasi-ld (quasi-2d) in a 

standard sense, i.e. both L^p = y/Dr^ and Lt = y/D/T exceed the thickness 
t and width w of the sample. However, the elastic mean free path I is smaller 
than t and w. Therefore the diffusion is three-dimensional at short times. 
Since the main contribution to the dephasing at T = comes from high 
frequencies, the 3d expression 
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applies. (The cut-off dependent numerical factor b is of order 1.) On the 
other hand, the temperature-dependent part of the rate, see Eq. (|l5|), comes 
from low frequencies or long times, when the diffusion is one-dimensional, 
and the corresponding form should be used. 

In Ref. ^the maximum dephasing times for ten (6 quasi-ld and 4 quasi- 
2d) samples fabricated from the same material (AuPd) with practically 
the same resistivity were found to be nearly universal t^°-^ (0.8 ... 2) x 
10~^^ sec, independent of the sample geometry. In the Saclay-Michigan 
experimentaHZH the maximum dephasing times were found to be 2 to 3 or- 
ders of magnitude longer than those in Ref. |9|. Furthermore, r^"^^ for the 
silver sample was about an order of magnitude higher than that for the cop- 
per sample in spite of similar (though not identical) parameterfH. Observing 
this difference the authors!^ suggested that the low temperature saturation 
of Tip in disordered metal wires is material dependent and not universal. 

For our comparison we will use the data for 8 samples of Ref. ^ with 
nominally identical resistivity p ~ 24 //O cm and D ~ 1.5 x 10~^ m^/sec 
(samples C to F and H to K), 2 samples of Ref. ^ {Ag and Cu) and one gold 
sample {Au) of Ref. |2^ (denoted there by AuMSU). Since for all samples of 
Ref. P the maximum dephasing times were nearly the same, we only quote an 
average r™''^. Experimental values of r™"^ and our theoretical predictions 
for T(^o (Eq. (^) with b = 1) are summarized in the following Table: 



sample 


r™'^^ (sec) 


T^o (sec) 


C to F and H to K (averaged)i 


1.3x10"^^ 


0.3x10"^^ 




10x10"'^ 


2x10-^ 




1.8x10"^ 


0.3x10"^ 


AuE 


8x10"^ 


4x10"^^ 



We note that for all the above samples the expression (35) with 6=1 gives 
correct order-of-magnitude estimates for r^°^, even though these times dif- 
fer by 2 to 3 orders of magnitude. By adjusting the cut-off parameter in the 
range b ~ 0.3 . . . 0.6 we could produce a quantitative agreement with all the 
mentioned experimental values. Alternatively we can eliminate the ambigu- 
ity related to the cutoff b by considering the ratios between the maximum 
dephasing times for different samples. One finds: 





experiment 


theory (eq. (|35|)) 


^AuPd /Ag 


1.3 X 10^3 


1.5 X 10^3 


vo / vo 


0.7 X 10"^ 


1 X 10"^ 


^AuFd I^Au 


1.6 X 10"^ 


0.8 X 10"^ 



This agreement strongly supports the conclusion about the universality 
of the low temperature saturation of t^^ in disordered conductors. To com- 
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plete this summary we mention that also our Id and 2d expressions for the 
maximum dephasing time ( p!^ ) give correct order-of-magnitude estimates for 
^max £qj, ^Yjq above samples because the mean free path is comparable to the 
width and thickness of the wires. 

The comparison with the standard theory! is much worse. The authors 
of Ref. ^ found that the value of r^p which follows from AAK for their 2d 
samples differs already at IK significantly from their experimental values. 
At the lowest temperature, T ~ 80 mK, this would imply a deviation by 
3 to 4 orders of magnitude. Only for the thinnest samples the observed 
temperature dependence was consistent with oc T^'^/^ . However, in these 
cases the values of the dephasing time were such that Tt^^ < 1 for the whole 
temperature range, whereas for the standard theoryS Tt<^ ^ 1 is required. 

We now analyze the temperature dependence of the dephasing time r<^ 
detected in the experiments of Refs. ^ and 27,28. For quasi-ld samples Eq. 
([TBI ) can be written in the following dimensionless formtJ 




where Tr 







^r%^V2D' 



(36) 



Here i?K = /i/e^ ~ 25.8 is the resistance quantum. For a comparison with 



6.0x10 



4.0x10 - 



2.0x10 - 




Fig. 2. Dephasing rate 1/t<^ measured in Ref. ^ for the samples A (squares), 
B (circles) and F (open symbols) as a function of temperature T. The 
straight lines are guides to the eye. 



the power-law temperature dependence predicted by AAK the experimental 
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data had been presented on a double- logarithmic plot. However, since the 
temperature Tq ~ lOK is high and the dependence of Eq. p^ ) is roughly 
linear, it is more useful to re-plot them on a linear scale. For samples with 
similar parameters p, D and t^q the slope for quasi-ld samples should depend 
only on the cross section tw. In particular, since w for the sample F of Ref. 
^ is reported to be 4 times larger than for the samples A to E, one expects 
the slope to differ by a factor 4. The experimental data for the samples A, B 
and F, presented in Fig. 1, are in good agreement with these expectations, 
with a ratio 4.05 between the slopes corresponding to the samples A and 
F. On the other hand, the magnitude of the slopes turns out to be larger 



than predicted by Eq. (35,^). The difference is as large as a factor of order 
8 when 6 = 1, reducing to a factor 3... 4 when the cut-off parameter is 
chosen in the range b = 0.3 . . . 0.5. This agreement is acceptable within the 
accuracy of our theory. We conclude that our theory accounts for the main 
observations of Ref. U. 



Now we will address the results of Refs. 27,B8. In Fig. 3 of Ref. 27 



the authors presented their data for Au, Ag and Cu samples. These three 
samples had similar classical parameters, but both the magnitude and the 
temperature dependence of r<^ differ drastically. This observation led the 
authors to conclude that the behavior of r^p and, in particular, its saturation 
at low temperature may be material dependent. Later it was demonstrated 
by the same grouped that the unusually low values of in goldE3 were most 
likely due to high concentration of magnetic impurities. A pure gold sample 
was fabricated and a similar behavior was observedc3 as had been found 
previouslyE3 for the Ag sample: showed no saturation on a log-log plot 
down to T ~ 50 mK. At the same time it was confirmecS that Cu samples 
with very high purity showed saturation similarly to the earlier Cu sampled. 

A material dependence of Tip cannot be ruled out in general. However, 
we find that such a dependence is not needed in order to quantitatively 
explain the seemingly different behavior of Ag (no saturation down to T ~ 50 
mK) and Cu (clear saturation already at T^TOO mK). The 'saturation' 
temperature Tq, defined in (^), sets the scale at which thermal and quantum 
contributions to the dephasing rate become comparable. It is obvious that 
for otherwise identical parameters the thinner samples will show saturation 
at lower temperatures. For the parameters of the Cm0, A(^ and Ai^ 
samples and with the cutoff parameter 5=1 one finds from Eq. (^^) 
T(f" PS 2.8 K, T^^ « 90 mK and Tg^" 40 mK. It is important to note 
that, for example, the ratio Tq'^/T^^ ~ 30 is not sensitive to ambiguities 
related to the cutoff parameter b. We conclude that the temperature where 
quantum effects set in can easily vary by more than an order of magnitude 
for samples with similar macroscopic parameters. 
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T, K 

Fig. 3. The values of t^{T) for the samples Cu (circles), Ag (open squares) 
of Ref. ^ and Au (black squares) of Ref. Dotted lines are the best fits 
to the formulal = AT^/^. Solid lines are the best fits to Eq. (|6|) with 
both Tipo and Tq varied. 



In Fig. ^ we replot the data points for three samples of Refs. ^J2^ on a 
linear scale. Best fits to the Eq. (^) are shown by solid lines, while dotted 
lines correspond to the standard theory!. Obviously the saturating Cn-data 
can be fitted only with a non-vanishing dephasing rate at T = 0. Also the 
'non-saturating' samples Ag and Au can be fitted with a finite t^Qj although 
one cannot draw definite conclusions since the temperatures Tq are too low. 
From the fits we find T^"" = 1.2 K, T^'^ = 14 mK, Tq^" = 16 mK and 
T^Q = 2 ns, = 30 ns, t^q = 20 ns. The deviations from the values given 
above can again be reduced if we adjust the cut-off parameter in the range 
b = 0.3 . . . 0.7. We conclude that the experimental data of Refs. 27,^ are 
consistent with the predicted low-temperature saturation of the dephasing. 
No material dependence of r<^ needs to be assumed to get this agreement. 
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